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Abstract 

We show that the free relativistic wave equation which describes the particle 
without or with rest mass has more than one part of energy spectrum. One 
part of energy spectrum is beginning with rest energy and it is not limited 
by above. This part of spectrum is called by us as normal. Another part of 
energy spectrum is beginning with the zero energy . This part of spectrum 
is called by us as anomalous since the zero energy corresponds to the infinite 
group velocity.The presence of the zero in the energy spectrum permits to 
consider the Bose-Einstein condensation. We show that the heat capasity has 
the finite discontinuity at the condensation temperature. The last means that 
we have the phase transition at the condensation point. 

I. I. INTRODUCTION 

It is well known that the construction of the thermodynamics requires the knowledge 
of particle energy spectrum because the principal thermodynamic characteristic , namely, 
equilibrium temperature distribution function depends on particle energy (frequency). We 
denote the last by f r , 

f r = (exp(huj — fi/hsT) + r) _1 , r = 0, ±1. 
where the relativistic energy Tluj is equal to 



hoj = hcJ k 2 + k"h 



Here k 2 is a square of 3-D momentum and ko is the inverse Compton wave length. We have 
superhigh energy at k » k Q , and in this case energy is equal to %uj ~ Tick (ultrarelativistic 
case). In the opposite case k « ko the energy spectrum of Eq(l) transforms into non 
relativistic energy spectrum. Energy spectrum of Eq(l) corresponds to eigenfunctions (plane 
waves) of Klein-Gordon -Fock (KGF) equation 

\I/ < - + ' ) = exp (z(— ut + k z z)) exp i(k\X\ + k 2 x 2 ). (2) 

Later we will explain the representation of the solution of Eq.(2) as a product of two 
exponential functions . The plane waves of Eq.(2) are the basis for Green functions (GF) 
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of KFG equation in the time-like Minkowski space M^ + \ It is known from long ago (see 
Stiickelberg, Feynman, Fierz [1-3]) that the construction of the causal relativistic propagator 
(CRP) requires the knowledge of GF in space-like Minkowski space as well. Considering 
relativistic quantum mechanics in so-called characteristic representation (CR) we can easily 
prove the last assertion [4]. Let us discuss the main principles of relativistic dynamics of 
CR in a brief form. In CR we have to single out one of the spatial coordinates in a wave 
operator (for the sake of definiteness we assume it to be z and let us write the relativistic 
wave equation in the form of 1-D telegraphic equation with the operator parameter a 2 , 

( d 2 Id 2 \ 

+ a 2 )^(x 1 ,x 2 ,z,t) = 0. (3) 



where 



dz 2 c 2 dt 2 



d 2 d 2 
a 2 = — + — - k 2 = A ± - hi 



dx\ dx 
is elliptic operator. 

Dealing with CR it is convenient to classify states by the eigenfunctions of the elliptic 
operator Aj_. This operator has two kinds of solutions. The first solution is a plane wave. 
The last gives the contribution to the GF in the time-like space M^ + \ The second solution 
is the Macdonald function K , 

A ± K (k ± x±) = k 2 ± K (k ± ,x±), 



where x± = \Jx 2 + x 2 and k± = \jk\ + k\. 

Therefore together with the plane wave there is one more solution 

¥~ ] = exp (i(-ut + k z z)) K (k±x±), (4) 

where uj is equal now to 



U = cy/k* -k 2 ± + k 2 . (5) 

We have shown in [5] that the both solutions of Eqs.(2) and (4) are the basis for two 
GF. One is defined inside the light cone and another is out of the light cone. It is especially 
simply to obtain two GF for photons when rest mass is equal to zero. Taking into account 
that the 1-D telegraphic equation (3) has the fundamental solution in the terms of Bessel 
function (BF) of the first kind of order zero we obtain [4] 

W-rf ,2W A ^ I n\ _ 1 ^1 ~ r / c ) _ 6 ^ ~ ^ 



Gff = </o (vW^Pj) |0> 



2tt cr 2tt 



where S(z) is Dirac delta function, and here it's taken into account that 5 (\t\ + r/c) — and 
also 



G# = Jo (y/(z*-M)(AS) 
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2tt (r 2 - c 2 t 2 ) 
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Here 



and 



|0> 



4tt 2 



expi(k 1 x 1 + k 2 x 2 ) = 8(x 1 )S(x 2 ) 



= 



2tt 



K (k±x±)k±dk± = 



2tt {x\ + x\) 



(6) 



are the initial values of GF which are given along the characteristics \t\ = \z/c\. The sum of 
these GF is CRP D c [2] 



2 ph 



ph ' 



The decomposition of the BF into the plane waves is presented in [5]. If rest mass is not 
equal to zero then the operator argument in the BF is changed: — Aj_ — > k% — A±. Inside the 
light cone we have the function = A s , where A s is the even solution of KFG equation 
[6] and outside the light cone we have 



G<-> = Jo (VC^-c^XA.L-fcg)) 



koK^koVr 2 - c 2 t 2 ) 



Vr 



2 t 2 



where 



with (6)) 



0) is the initial value of propagator G^ ) along characteristics \t\ = \z/c\ (compare 



2tt 



0) = — / K (k±x±)k±dk_i 



ko 



k KijkoXj) 

2/T Xj_ 



(7) 



Let us write the group velocities along the determinate direction for energy spectra of 
Eqs. (1) and (5). For Eq.(l) 



duj 



ck-, 



dk z ^2 + q 2 ± + k 2 



for Eq.(5) 



duj 



ck? 



dk * yjk* -k 2 ± + k 2 



> \c\ . 



What differs the spectrum of Eq.(5) from the spectrum of Eq.(l)? The main difference lies 
in the presence of minimal zeroes energy for which the velocity i^is maximal and goes to 
infinity. The zero energy (frequency) is realized if the transverse momentum k± ^ and is 
equal to 



k± = \J k\ + k\ > ko, at uj — 0. 
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Due to such extraordinary property that the infinitely- moving particle has zero energy, the 
spectrum of Eq. (5) is called by us as anomalous, although the definition of this spectrum as 
anomalous is not good. The state describing this spectrum is needed as well as normal. As 
it has been mentioned above, the CRP -D c consists of these two states (see Eqs.(2)and (4)). 

Energy spectrum of Eqs.(l) and (5) (dispersion curves) are mapped on the plane (w/c, k z ) 
by half-hyperbola which belong to different quadrants of the plane (u/c,k z ) divided by 
asymptotes \ou/c\ = \k z \ . 

In this paper we examine the special features of thermodynamic characteristics concerning 
statistical integral, internal and free energies, enthropy etc. referring to relativistic gas in 
the anomalous state. 



II. THERMODYNAMIC CHARACTERISTICS OF THE NORMAL STATES OF 

RELATIVISTIC GAS 

For comparison we present in a brief form the known results of relativistic thermodynam- 
ics (RT) based on the spectrum of Eq.(l) [ 7 ] .First we consider the Maxwell -Boltzmann(M-B) 
statistics , that is r = .We denote the particle statistic integral in M-B statistics by Z 

Z = J rf 3 A;exp = ^klJeV^le W (^) de, (8) 

where = mc 2 jks is typical parameter of RT having the temperature dimensionality and 
ks is Boltzmann constant. If m is electron mass then 6 10 10 ° K. 

Let us denote the factor before the exponential function in integral of Eq. (8) by g(e) 

g(e) = 4nk 3 eVe 2 - 1, 

or in dimensional units (e = Tiuj/mc 2 ) 

Tiuj 



g(u/c) = 4nk { 



3 

mc 2 



mc 2 



~ I- (9) 



Eq.(9) gives the state density in the energy scale which is determined up to the factor 
g^3, where VL is the volume of relativistic gas. This density is equal to zero if Tiuj < mc 2 
or in temperature units T < ©.Due to such temperature restriction it is said that RT 
has academic interest and it can be applicable in astrophysical phenomena where there 
are superhigh temperatures and densities. The integral of Eq.(8) is easily calculated by 

oo 

substitution e = coshr and with the help of the table integral / dz exp(—q cosh z) = K (q). 

o 

It is respectively equal to 

Z = ^ ^S1 (10 ) 



where Ki is the Macdonald function of the second kind 

mTnordrnrD urnon 

T 



Let us calculate the asymptotics of Eq.(lO) at low temperature when !p >> 1 
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Z^(2vr) 3 / 2 A; 3 (|)lexp 



Statistic integral of Eq.(lO) does not depend on mass at high temperature when (^) << 1 



Let us compute the internal energy density 



T > 

3 



U = (hu} = J cPkhujexp ^ 



-hu>\ = kBT2 dZ_ 



k B T J dT 

It is convenient to find the internal energy in per one particle so-called middle energy 

T (hu) t2 d 2 T K l + K^ 

U = — = k B T-\nZ = »c + 

At low temperature the middle energy and heat capacity (HC) are equal to 



U ^/ic 2 + 3/2k B T, c v = %k B . 
The middle energy and HC at high temperature (®) << 1 are equal to 

U ~ 3k B T, c v = 3k B . 

If the relativistic gas obeys to Bose-Einstein (B-S) statistics that is r = — 1 then all thermo- 
dynamic characteristics can be found from the distribution 

JV=4rt °S — wjt) — ■ 

III. THERMODYNAMICS OF ANOMALOUS STATES 

Let us find the state density in an energy scale for the energy spectrum of Eq.(5). It is 
convenient to introduce dimensionless variables, namely, dimensionless energy (frequency) 
and dimensionless momentum variables 

e = u/k c, q z = k z /k , q± = k±/k . 

We write dispersion equation (5) in dimensionless variables 



In fact, it is necessary to make the transition from the momentum representation to energy 
representation with the help of transformation formulae which depend on the connection 
between full energy (hcu ) and rest energy (mc 2 ). It means that in dimensionless variables 
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the value e can be more than one or less than one. Let e > 1 then transformation formulae 
have the form 

q z = \J e 2 — 1 cosh v, q\ = y/e 2 — 1 cos u sinh v, qi = sin w sinh -y. 

The infinitely -small volume element d 3 k (d 3 q in dimensionless variables ) in the new variables 
has the form 

d 3 q = \I\ dedudv, (11) 

where I is Jacobian of transformation . It is equal to / = e\J e 2 — 1 sinht>. 

To obtain the state density in energy scale it is needed to integrate Eq.(ll) over two 
angles u and v .The integration over the angle u gives the factor 2tt. To integrate Eq.(ll) 
over hyperbolic angle v it is necessary to know its limits of integration. They are found from 
formula for the domain of the definition of transverse momentum k± ( q±in dimensionless 
variables ) 



Q± = ^ll + ll = Ve 2 - lsinhv. (12) 

The low limit of the hyperbolic angle v corresponds to the minimal value q±. Since 
g mm = i( r in dimension variables A;™ 111 = k (see Eq.(7)) then 



sinhw min = ,cosh^ min = \Jl + sinh 



2 v„ 



\]E l — 1 V£ — 1 

If the upper limit of the value q± is equal to infinity then the upper limit of hyperbolic angle v 
is equal to infinity too. In this case Eq.(ll) becomes infinite after integration over hyperbolic 
angle v .Therefore we cut the momentum k± by the value /c™ ax . It does not influence on the 
state with zero energy because 



u = at k™ x = ^J(kf^) 2 + kl > k or g max > 1. 
We can find the upper integral limit from Eq.(12) 

qT x u ^J(qT x ) 2 + e 2 -l 

smhv max = cosh t; max = . 

y/e 2 — 1 \Je 2 — 1 

Integrating Eq.(ll) over the angles u and v we obtain the state density in the energy scale 

g(e)de = 2nk% (e^e 2 + ( g max )2 _ i _ e 2 ^ de (13) 

Further we consider the case e < 1. It means that in the dimensional variables energy is 
less than rest energy mc 2 (or e < 1). In this case the transformation formulae have the form 



q z = \/\ — e 2 sinh v , q\ = Vl — £ 2 cos u cosh v , 
q2 = \J\ — e 2 sin u cosh v , 

and the volume element in dimensionless variables has the form 
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d 2 q = \I\ dedudv = ey/l — e 2 dedud(smhv) (14) 



where / = eyl — e 2 cosh v is Jacobian of transformation. 
We can find the limits of the angle v from the expression 



Q± — \Jqi + Q2 — Vl — s 2 coshf . 
Repeating preceding reasonings we find 



cosh v min = , 1 = , sinh v m ™ = 



VT^' mm VT^ 

and 



cosh v mSLX = j= , sinh v 



• /(g™ ax ) 2 + e 2 - 1 



, . max . 

Performing in Eq.(14) the integration over angle variables we find the state density coinciding 
with Eq.(13). 

Let us rewrite Eq.(13) in dimensional variables 



If the value k in Eq.(15) goes to zero then the form of this equation is not changed. 
From Eq.(15) it follows 



\img(uj/c) = nko ((k™ x ) 2 - fcg) at u -> oo. 



Introducing the notion of effective mass m* = my(g™ ax ) 2 — 1, we rewrite the state density 
of Eq.(13) in a more compact form 



«U O \ 2 • (— V s 2 | (16) 



Now we can consider the statistic integral of anomalous state in M-B statistics in energy 
representation 



Z = 2nkl j Us 2 +{^) 2 -eA exp (-^) de (17) 



One integral in Eq.(17) is calculated in an elementary manner 



oo 



e 2 exp ( -^f) de = l( 1 - 



T J V0 
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The second integral of Eq.(17) is calculated with the help of substitution e = ^- sinh 6 and of 

oo 

table integral / d6 exp(— q sinh 6) = jA (q), where A Q = H — N and H is Struve function 

o ' " ' 

and iVo is the Bessel function of the second kind of the order zero. From what has been said 

above it follows that 



* 2 2 I 

where 6* = = J (g™ ax ) 2 _ \ j s effective temperature and 

5(g) = 2^/g 2 - A /q - 4/irq 3 , A, = H, - TV-, i = 1, 2 (19) 

The spectrum of Eq.(5) unlike the spectrum of Eq.(l) permits to consider B-E condensation 
of relativistic gas. In this case the cutting parameter /c™ ax (or effective mass to*) can be 
expressed in the terms of condensation temperature T c and gas particle number N a per unit 
volume. The condensation temperature T c is defined from the condition of vanishing of 
chemical potential /x in the gas particle density N a . From the point of view of B-E statistics 
the N a is equal to 

= £ B (f!) (20) 

where k£ = rn jf L and the value B(q) is defined from Eq.(19). At low temperature the 
condition >> 1 is always performed and we can use the asymptotics of the function 
B(q) w 2/nq 2 at q » 1 , 

iV a « 2n,(2)(k* f = 2tt ? (2)A$ (21) (21) 

where q(s) is zeta-function of Riemann . 

From Eq.(21) it follows that the effective mass to* satisfies the condition 

to*T c = const. 

The gas particle density iV with the zero energy (condensation particles) and the density of 
overcondesation particles Aqof relativistic gas are equal to 



The internal energy density U of relativistic gas in temperature interval < T < T c is 
defined by integral 

t, = ^ / = Ajvv £ D (rfl) (22) 

S eX P J ~ 1 "=1 V 7 

where = 2/7rg - 12/vrg 4 - Ai/g - 3A /g 2 + e^/g 3 . 
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Taking into account that at low temperature (q » 1) the asymptotics of D(q)— function 
is~ 4/7rg 3 we have from Eq.(22) 

U » 4^(3)A;* 3 m*c 2 (J^f = l,5N a k B T (^j (23) 

In Eq. (23) the approximate equality has been used 2<j(3)/<t(2) = 1,5. 

The HC C v in the temperature interval < T < T c is equal respectively to 



C v = A,5N a k B (-) 



Let us investigate the relativistic thermodynamics in the temperature interval T — T c « 
T c . We take into account that N is equal to zero and \fx\ /k B T « 1. Then the particle 
density in relativistic gas has the form 

n ° - 7 exp( 4^)-i k M *°* )M2) (^ exp © (24) 

In Eq.(24) we have used the asymptotics of B(q) and the approximate summation with 
the help of Euler formula £ exp ^" g) « <j(2) exp(-g) which is valid at q « 1. 

n=l 

Substituting Eq.(21) in the left side of Eq.(24) we obtain the equation for the chemical 
potential \i , 

(T\ 2 , ,, _ N _ ( T ' 



\T J 

\J- c / 



exp(^/k B T) -> /i = 2/c B Tln ^— j at T > T c 



All the derivatives of chemical potential /x are discontinuous at the condensation temper- 
ature T = T c , in contrary to non relativistic thermodynamics where the first derivative of 
chemical potential \i is continuous at the condensation temperature T = T c , 

/i'(T c -0) = , /i'(T c + 0) = 2k B , 

and etc. 

Let us calculate the internal energy density U in the temperature interval T > T c , 

oo 



U = mc 2 J , 

)l T 



exp 



47rA;* 3 m*c\(3) (£) exp (-^) ~ 1, 5N a k B T, 

and HC is equal to C v = l,5N a k B . Thus, in the temperature interval T > T c the internal 
energy and HC of relativistic gas in the anomalous state coincide with energy and HC of 
ideal non relativistic gas. 

Let us write the expression for full internal energy using the Heaviside step-function 9{t) 

i 

U{T) = l,hN a k B T\^06{T c -T) + 6{T-T c )} , 
U(T C ) = 1, 5N a k B T c . 
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At the condensation temperature T = T c we have the heat capacity jump 



C v (T c -0)-C v (T e + 0) _ 

C v (T c + 0) " (25) 

Eq.(25) does not contain the arbitrary parameters and therefore it admits the experimental 
verification. 

Enthropy variation can be found from the expression dS = whence , after integration 
we obtain 

S{T) = h B N a | Vc " T) + (in (|) ^ + l) 0(T - T c ) } 
and S(T C ) = \N a k B . 

We can find the expressions for free energy F = U — TS, for an enthalpy H = U + PV = 
TS + fiN a 9(T — T c ), the equation of state of relativistic gas etc. These main values have the 
fracture (or angle point) at the condensation temperature T c . It means that the derivative of 
these values with respect to temperature T has the discontinuity of the first kind at T = T c . 
Therefore these thermodynamic characteristics cannot be expanded into Taylor series in the 
vicinity of the critical point T = T c . 

For comparison let us write the particle density N ntT of non relativistic gas at the con- 
densation temperature T c , 

,rp . 3/2 

N n , r = {2nfl\{mkl (^) 
The ratio of two densities (see Eq.(21)) 

N a ?(2) m* (T c \ 1 ' 2 



N n>r v^(3/2) m 19 

shows that at low temperature T c the numbers N a « N n>r . This inequality strengthens 
if m* << m. 



IV. CONCLUSION 

Our reasonings are at least in harmony with the papers of Bilaniuk, Terletskiy , Feinberg 
(see ref.[8] and refs. therein ) in which it is proved that the relativistic quantum mechanics 
is incomplete if we are restricted by energy spectrum of Eq.(l). We know that the RCP 
D c consists of the solutions of different nature (oscillatory and exponential damping). The 
last corresponds to energy spectrum of Eq.(5) which we call anomalous. We could call this 
spectrum as tachyon although tachyon term is referred to the solutions of KGF equation 
with imaginary mass[9]. We believe that the energy spectrum of Eq.(5) can be applied 
to the description of thermodynamics of models with the observed phase transition of the 
second kind. 

REFERENCES 



10 



1 E. Stiickelberg, Helv.Phys.Acta, 14, 588, (1941) 

2 R.Feynman, Phys.Rev., 76, 769 (1949) 

3 M.Fierz , Helv.Phys.Acta, 23, 731 (1950) 

4 A. Borghardt and D. Karpenko , J. Nonlinear Math. Phys., 5, 357 (1998) 

5 A. Borghardt, D. Karpenko, JETP, 85(4), 635, (1997) 

6 N.Bogoliubov and D. Shirkov, Introduction to the Theory of Quantized Fields (Inter- 
science, NY,1959) 

7 W.Pauli, Relativitdtstheorie, Ens. der Math. Wiss., 5, 539 (1921) 

8 O.Bilaniuk, V.Deshpande, E.Sudarshan, Amer. J.Phys., 30, 718,(1962) 

9 G.Feinberg, Phys. Rev., 159, 1089 (1967) 



11 



